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Abstract 



We simplify and generalize an approach proposed by Di Vecchia and Ravn- 
dal to describe a massive Dirac particle in external vector and scalar fields. 
Two different path integral representations for the propagator are derived 
systematically without the usual five-dimensional extension and shown to be 
equivalent due to the supersymmetry of the action. They correspond to a 
projection on the mass of the particle either continuously or at the end of 
the time evolution. It is shown that the supersymmetry transformations are 
generated by shifting and scaling the supertimes and the invariant difference 
of two supertimes is given for the general case. A nonrelativistic reduction of 
the relativistic propagator leads to a three-dimensional path integral with the 
usual Pauli Hamiltonian. By integrating out the photons we obtain the effec- 
tive action for quenched QED and use it to derive the gauge-transformation 
properties of the general Green function of the theory. 
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I. INTRODUCTION 



The problem how to describe spin in a path integral has a long and twisted history. 
This is mostly due to the fact that a path integral is determined by the classical Lagrangian 
(or Hamiltonian) and a classical analogue for the internal spin of a particle is not readily 
available. Martin [|lj apparently first suggested to use anticommuting Grassmann variables 
for this purpose. This can be made plausible when one recalls that the spin operator s = 
ha/ 2 of an electron fulfills 

in the classical limit. Consequently one can describe a spinning particle by its bosonic part, 
the usual trajectory x(t), and a fermionic degree of freedom given by a Grassmann valued 
function ((t) Q Brink et al. Q noted an important supersymmetry between the bosonic 
and fermionic parts of a relativistic massless Dirac particle and Berezin and Marinov ||] 
showed that massive particles can be described by adding a fifth component, (sit), to the 
spin variable. The reason for this peculiar addition is that in the rest frame of the particle 
the spin is intrinsically 3-dimensional (cf. Eq.(|l|)) and a covariant 4-dimensional description 
therefore has superfluous degrees of freedom which have to be cancelled by the fifth spin 
variable ||. 

There is now a vast literature about spin in path integrals (a partial list of references is 
[6-10]) which discusses various aspects of this approach. In particular, Fradkin and Gitman 
HTT|] have given a straightforward way of constructing the corresponding relativistic propaga- 
tor. In addition to the dependence on bosonic and fermionic trajectories mentioned above, 
their formulation has the special feature that as well as the usual Schwinger proper time 
a Grassmannian partner to it is required. Representing Dirac particles in a first quantized 
form in the "worldline formalism" has become popular for perturbative calculations in QED 
and QCD [12-14]. These one-loop calculations of the effective action are simplified by the 
fact that only Green functions on a circle (with simpler boundary conditions) are needed. 
More recently, the method has also been used in order to derive derivative expansions of the 
one-loop effective action in 2 + 1 and 3 + 1 dimensional QED Ji~5|j . 

Although sufficient for many purposes the Berezin-Marinov introduction of the fifth 
spin variable is an awkward one: there is no clear physical picture associated with it and 



the corresponding multiplication of the propagator with the Dirac matrix 75 [1.1 fl is very 



unnatural in a parity conserving theory. A four-dimensional approach, which did not get 
very much attention up to now, is that proposed by Di Vecchia and Ravndal in which 

the unwanted spin degrees of freedom are simply projected out Q. 



x It should be noted that there are other approaches, e.g. using coherent state path integrals [|j, 
which we will not consider here. 

2 During the course of this work there appeared a publication |[^] in which the fifth spin variable 
is also eliminated but by a different nonlinear technique. 
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It is the purpose of the present paper to develop this latter approach further and to show 
that it has attractive features. In particular, in Section 2, we will calculate the propagator 
for a Dirac particle in an external vector field and demonstrate that the projection mentioned 
above can be done in two different ways: either at each time step during the evolution of 
the system or at the end. We will refer to the former as the "local" projection method and 
to the latter as the "global" projection method. In Section 3 both procedures are shown to 
be equivalent due to the supersymmetry between bosonic and fermionic variables. However, 
the global projection leads, in general, to simpler expressions without a Grassmann proper 
time. In addition, an inherent coupling between orbital and spin parts which is already 
present for a free particle is removed by the global projection method. Section 4 contains 
the nonrelativistic reduction where we start directly from the path integral representation of 
the Dirac propagator and show that this reduces to the three-dimensional spin path integral 
of the nonrelativistic theory. This is to be contrasted with Ref. || where the nonrelativistic 
propagator was derived starting with the nonrelativistic Hamiltonian and introducing three- 
dimensional Grassmann variables instead of obtaining it from the relativistic path integral for 
the Dirac propagator. In Section 5 we show that one can also describe a Yukawa interaction 
of the fermion (i.e. the particle in an external scalar field) in such a four- dimensional 
framework. As an application we derive the effective action in Quantum Electrodynamics 
in Section 6 and finally we summarize our results. 

Since we aim in making this paper self-contained we include an Appendix with a deriva- 
tion of the spin path integral which is somewhat different, more explicit and simpler than 
the one given by Fradkin and Gitman. Our conventions follow Bjorken and Drell |19| and 
in general we use the letters a — l and £ — ui to denote Grassmann variables, with some 
exceptions to comply with standard notation found in the literature. 

II. DIRAC PROPAGATOR IN AN EXTERNAL VECTOR POTENTIAL 

We are looking for the path integral representation for the propagator of a Dirac particle 

1 



G(x,y) = (x 



- gjjL(x) - M + i0 



y) (2) 



in an external field A^(x) where throughout this paper quantum mechanical operators are 
denoted by "hats" over the corresponding symbols. In the spinless (bosonic) case this can 
be achieved by using Schwinger's proper time representation for the quantum mechanical 
resolvent 



1 . f°° 



E - H + iO 



/ ciTexp i(E-H + iO)T . (3) 
Jo 1 J 



However, in the fermionic case we have to make sure that the operator H which plays the 
role of a Hamiltonian for the quantum mechanical system contains an even number of Dirac 
matrices. This is because in the classical limit (and also in the path integral) only an even, 
commuting object can represent a physical quantity. 



Fradkin and Gitman |TT[ achieved this by multiplying numerator and denominator in 



Eq. (g) by 75 and extending the Dirac algebra to five dimensions. However, it is much 



simpler to use the representation of Di Vecchia and Ravndal where we write 
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fl+ M j — 

fl — M 2 + iO 



fl-M + iO 

It is now possible either to exponentiate only the denominator which gives 



1 



dT (Jfl+ M) exp 



fl-M + iO 2k Jo 
or both numerator and denominator leading to 



iM 2 T \ / i ±3 
2k J 6XP V 



1 



fl-M + iO J ° 



dT / dx exp 



2k 



(M 2 T + M X ) 



exp 



2k 



The latter only holds if fl commutes with fl which is proved in Ref. [IT). Here 

ft" = - gtt\x) , 
and the Berezin integrals over Grassmann variables are defined as [EH] 



/ dx = , J dx X = !• 



(4) 



(5) 



(6) 



(7) 



(8) 



k,q is a parameter which reparametrizes the proper times T — > k T , x — > k X without 
changing the physics and is a remnant of the local reparametrization invariance of the action. 
It is thus related to the "einbein" |J. x is either called a (one-dimensional) "gravitino" field 
or, more appropriately in the present context, as the supersymmetric partner of the proper 
time, the "supertime". 

The Di Vecchia-Ravndal representation has several advantages compared to the stan- 
dard Berezin-Marinov form Q for the description of a massive spinning particle: no five- 
dimensional extension and multiplication with 75 are necessary and, as we will see in Sec- 
tion 3, the supersymmetric transformations are much simpler and more transparent. It can 
be considered as the result obtained by integrating out the fifth spin variable. A certain 
disadvantage is that not all exponents in Eq. @ are Grassmann even. The odd term 

exp r 2^ x ) (9) 

is to be considered as part of an operator which projects out fl— M fl7j and not as part of 
the evolution operator. In Eq. (|5[) this projection is done at the end ("global") whereas in 
Eq. (H) it is done at each time step during the evolution ("local"). 

It is essential that in both procedures the "Hamiltonian" which governs the proper-time 
evolution is even. In the global projection method it is given by 



w(n,x, 7 ) 



1 $ 
-fi 



n 2 



2kq 2kq 4ko 

whereas for the local projection method it reads 

W(n,£,7) = H(fl,x,j) 



d v A^ 



2k T 



ti.rx- 



(10) 



(11) 



In both cases the parameter kq can be interpreted as the "mass" of the quantum-mechanical 
particle. 
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A. Global projection 



We will first consider the projection after the end of the evolution, i.e 

iM 2 T^ 



G(x,y) 



2k 

i 

2«o 



d z (x 



i> - g4ix) + M 
( i$ x - g4{x) + M ) J dT exp 



dT exp y 

( iM 2 T 
{ 2^T 



2k 



exp 



( -iHT 



x 



exp ( —iH T 



V 

12) 



The remaining proper-time evolution operator can be written in path integral form following 
Fradkin and Gitman [11] but staying within a four- dimensional framework. In Appendix A 
we show that 



x 



exp ( — iH T 



Q \ rx(T)=x 



a) = ex P y ■ — 
• exp (Vo 



. , VxVpVi N sp ' m 

dT J Jx(0)=y 
T 

dt - v ■ x - n{ n, x, 2£ + r) 



(13) 



r=o 



where A spm is a normalization factor for the four- dimensional spin path integral as given in 
Eq. (A19) and we use antiperiodic boundary conditions for the spin variable £(t) 



Eq. ( |T3"D can be further simplified by shifting to the new spin variables 
so that the boundary condition becomes 

C M (o) + UT) = r„. 



(14) 



(15) 



(16) 



This introduces an additional boundary term: — |r • [C(T) — C(0) ] = C(T) ' C(0) • After 
shifting to the momentum (]?]) as integration variable we obtain 

G(x, y) = - 94(x) + M) exp ( 7 • — j J dT exp (-^-M 2 t) A spin J VxVILV{ 



exp C(0)-CCO +i dt 



i( ■ C - (n + gA{x)) ■ x -H{IL,x,2( 



• (17) 



r=o 



As usual we can perform the functional Il-integration since the Hamiltonian is at most 
quadratic in the kinematical momentum. We then obtain the final expression 



G(x,y) 



2nr 



( ife - g4(x) + M ) exp U • A j j™ dT N(T) exp ( - J-M 2 T 

VxV( exp{iS[x,(] } r=Q 



(18) 
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where 



N(T) 



V( exp C(0) • CCO 



dt(-( 



VU exp i 



dt 



2k 



(19) 



provides the proper normalization and 



S[x,C] = 



T 



dtL(x,x,(,Q - K(P)-C(T) 
^x 2 +i(.{- gx . A{x) - ^F, u (x) (»C 



(20) 



are the action and the Lagrangian, respectively. The first two terms in Eq. (p0|) correspond, 
respectively, to contributions from the orbital and spin degrees of freedom to the kinetic 
energy, while the last two terms are the contributions of the photon field coupling to both 
the electron's convection current and its spin current. The canonically conjugate momenta 
are given by 



and 



dL 

dx^ 



-k x^ - gA^x) 



dL 



so that the canonical Hamiltonian becomes 



H 



dL 



ig 



(21) 



(22) 



(23) 



when expressed in terms of (generalized) coordinates and velocities. In terms of coordinates 
and momenta we have the relation H = TC\ p ^_ p . This is a consequence of our metric which 
gives exp(— ip ■ x) as plane wave and therefore leads to the form / dt [— p ■ x — H] for the 
action in the phase space path integral (jT3|). 

The free Dirac propagator in momentum space is readily obtained from fll8|) since orbital 
and spin variables decouple. The £-path integral cancels against the normalization factor 
jyspm an( j x-path integral gives just the usual free bosonic evolution kernel. Thus 



G<°>(p) 



d 4 x e ip - x 
d A k 



2kq 



+ M) I dT exp 



2kq 



-M T 



(27T) 



e ifc ^exp 



2k 



-k 2 T 



p* 



M 2 + i0 p'-M + iO 



(24) 



independent of the reparametrization parameter kq. 
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B. Local projection 



The local projection method follows along the same lines with two differences: first we 
have an additional integration over the supertime x an d second there is an extra term in 
the action due to the additional term in Eq. ([□]). Thus 



G \x, y) = exp (7 • ^ ) dT N(T) J d x exp 



1 -r2 



(M T + M X ) 



2k,q 

JvxV(exv{iS'[x,(}} T=0 (25) 



with 



S'[x,C}= f dtL'(x,x,(,(; X ) - i((0)-((T) 
Jo 



L '(x, x, c, C; x) = L{x, x, C, + y X (26) 
= - ^x 2 + i( ■ C + ^ x ■ CX - gi ■ A{x) - ^F^x) ■ 

Note that there is now a coupling between orbital movement and spin, even for the free 
particle. This is the same mechanism which at high energy aligns the spin of a Dirac 
particle along (or opposite) to the momentum whereas a non-relativistic particle with spin 
is unaffected. 

Since the spin degrees of freedom appear at most quadratically it is also possible to 
integrate them out completely and reduce the path integral to a bosonic one modified by 
a "spin factor" |2T|j2^|. The price to be paid is that this spin factor is highly nonlinear in 
the external fields. This prevents an analytic integration over the boson fields to obtain an 
effective interaction for the fermion only, as is done in Section 6. 



III. BOSONIC AND FERMIONIC TRANSFORMATIONS 

We next discuss the transformation properties of the Lagrange function in the local 
formulation |3,23|]. The corresponding ones for the global formulation can be obtained by 



setting x = . There are two kinds of transformations which leave the Lagrange function 
L' in Eq. (|2~6D invariant (up to a total derivative): 

(i) bosonic transformations (reparametrizations) 

= b(t) 

6(» = b(t) C (27) 

2 d ( b{t) \ 



°dt 



«0 



6L'= j t [b(t)L>] , (28) 
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where b(t) is the infinitesimal parameter of the transformation which, in principle, 
could have an arbitrary time- dependence (local transformations). However, for quan- 
tization the reparametrization "gauge" has to be fixed 0,0] which in our case, by 
construction, was taken to be k = constant. This means that we only can allow 
5k = constant, or 



b(t) = bo + ht 



(29) 



Note that b\ = corresponds to proper time translations, e.g. 5x^ = x fl (t+bo)—x fl (t) 
&o£ M + ..., and 6o = to proper time scalings, e.g. Sx^ = x^{t+bit)—x^{t) = bitx^ 1 + . 

(ii) fermionic (supersymmetric) transformations 



6x» = ia(t) C p 
*C = f a(t) 



x' L 



K T 



Cx 



5k 

$X 
5L' 



T 



a(t)x 



-iKoTa(t) 

|[a(t)(-|*.C + M. C 



(30) 
(31) 



where a(t) is the infinitesimal Grassmannian parameter of the transformation. Again, 
since k,q and x are by construction time-independent, one can only allow transforma- 
tions with a(t) = constant. For x — Ravndal [|1^] has shown that similar to the 
bosonic case it is also possible to generate the fermionic transformations by a shift in 
the proper time if a Grassmannian partner of the proper time t is added. This allows 
for a concise supersymmetric formulation of the action. In this Section we will show 
the generalization of Ravndal's transformations to the case with spin-orbit coupling, 
which includes a special scaling of the "supertime" in addition to a shift. 



Since the change of the full Lagrange function is a total derivative, Noether's theorem |2~3 



allows us to define quantities which are conserved classically. For the bosonic transformation 
with bi = (i.e. proper time translations) we have Skq = and Eq. ( p8|) therefore leads 
to the conservation of the canonical Hamiltonian (p3|). It can be shown that proper time 
scalings (po = 0) lead to the same result 0. For the fermionic transformations we find from 
Eq. (|3lD that the projection of the spin variable on the kinematical momentum — n x 



Q 



-KqX^C 



(32) 



is conserved classically without the spin-orbit term 17|. 



3 In this case the reparametrization parameter kq, which is also changed, is not a dynamical 
variable for which the equations of motion can be used. Consequently the change of the corre- 
sponding Noether charge Q = tH with time is proportional to K^dL/dntQ = H, which gives no new 
information. 



S 



Quantum mechanically the Noether charges either become conserved operators or, in the 
functional formalism, their conservation implies that certain averages, i.e. Green functions 
with the Noether charges as insertions, stay time- independent. For quantum mechanical 
averages we will use the following notation 



(0) s = exp 7 



dT 



Vx V(<D(x,() e 



iS[x,C] 



r=o 



(33) 



To be specific, we consider the fermionic transformations with \ = because their Noether 
charge fl3~2|) does not dependent (explicitly) on the interaction and we make a local, time- 
dependent transformation 



x(t) = x'(t)+ia(t) C(t) 



C{t) = C(t) + ^ a(t) x(t) 



(34) 



in the path integral. We assume that a(0) = u(T) = so that we do not have to consider 
boundary contributions. The Jacobian for this transformation is 1 + 0(a 2 ). Since the path 
integral does not change its value we obtain (omitting the primes) 







dt 5L 



(35) 



where 



5L = % a(t) 



d 
dt 



-^x-C-gA-C } +i<\(t) 



(36) 



The first term is what we obtain for a global, time-independent transformation in Eq. (pip. 
Performing an integration by parts in the second term (no boundary terms) the result is 
then 







dt a(t) 



(37) 



or since a(t) is arbitrary 



j t ( -K x{t) ■ ((t) j ^ 



(38) 



for all times. 



A. Supersymmetric formulation 

It is convenient to write the Lagrange function for a relativistic spinning particle in 
explicit supersymmetric form by combining orbital and spin degrees of freedom into a "su- 
perfield" 0] or "superposition" [[HJ 

X»(t,9) = x"(t) + a6(»(t) . (39) 



9 



Here 9 is an additional time-independent Grassmann variable which acts as a superpartner 
of the proper time t and a a suitably chosen constant. If, in addition, a "superderivative" is 
defined as 



D 



d__ d_ 
89 dt 



then 



(W 



K 



DX,D 2 X» = -^x 2 + i(-( 



(40) 



(41) 



generates all terms in the free Lagrangian of the spinning particle provided the constant a 
is chosen as 



(42) 



Note that in this compact form only first-order derivatives appear since D 2 = —d/dt. 
In the local projection approach one needs an additional x-dependent factor P,|13| 



e(9 X ) = 1 + -^9 X 



(43) 



in the integrand of Eq. (|41|) to account for the explicit spin-orbit coupling. Thus the corre- 
sponding free action is 



S '[X] 



dt / dd e(0x) 



2 



DX ■ D 2 X 



dt 



— w x + < ■ C + tttC ■ xx 



1 

T~ 



(44) 



The interaction of the Dirac particle with an electromagnetic field takes the equally simple 
form 



L e . m . = 9 J d9DX^{X) 



(45) 



which is easily proved by expanding the "superposition" X and performing the Berezin inte- 
gration. Eq. (flop thus contains both the convection current and the spin current interaction. 

For x = Ravndal and Di Vecchia P^JTTH have given a simple way of generating both 
the bosonic (with bi = 0) as well as the fermionic transformations by a shift in the proper 
times t and 6: 



t -> t + b + e 9 
9^9 + e , 

where e and bo are constants which may be zero. Indeed, the superfield changes into 



(46) 



X(t, 9) -> X '(t, 9) = x(t + b + e0) + a(9 + e)((t + b + et 



x + b x + ae( + a9 



C + b ( - -x 
a 



+ 



(47) 
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and if we set 



e = -a (48) 
a 

we obtain both transformations ( p7|) and (|30|) for the individual components of the superfield 
in the special case x = • This is not only more transparent but also treats bosonic and 
fermionic transformations on an equal footing. The equations of motion and the conserved 
quantities can also be formulated compactly in this formalism. 

We can generalize the transformations ([16]) to the case \ 7^ by observing that any change 
in t, 9 leaves k , \ unchanged, since these quantities are by construction time-independent. 
This means that necessarily 

5k = (49) 
6 X = 0. (50) 

While the latter condition is fulfilled by a constant parameter a in the fermionic transfor- 
mation (see Eq. (|30"D ) the former one requires that the bosonic scaling parameter bi is not 
arbitrary but given by 

h = ^ax- (51) 

Using Eq. ([48]) we then find that 



t'= (l + ^ex)t + b + e6 

generate the x-dependent supersymmetric transformations with 5k = 5\ = . Although 
this constitutes a scaling of the bosonic time t by a factor 

l=(l + ^x) (53) 

the fermionic time 9 is only scaled by yl . Consequently D scales by 1/ \fl . Since the 
spin-orbit factor ([43]) scales again with I and the Berezin integral over 9 transforms inversely 
compared to a bosonic one the free action is easily found to be invariant under scaling. 
We also note that for two times t\, £2, 0i, #2 the combination 

T 12 = h ~ t2 + 9 X 9 2 (54) 
y/e(9 lX ) e(9 2X ) 

is invariant under the shift and scaling (|5^) of proper times. This is the generalization of 
a result which is well known for \ = [ETJ and is important for extensions of the polaron 



variational approach to QED [26 
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B. Equivalence of local and global projection 



We are now able to prove the equivalence between the local projection method and the 
global one. We give here a somewhat different and more explicit derivation than the one 
sketched in Ref. [27|. We start from the local formulation and perform the x ~ integration. 
This gives 



G'(x,y) = exp (7 • ^ 



x(T)=x 
x(0)=y 



Vx 



dT N{T) exp 
V( e iS 



C(o)+c(T)=r 



2kq 

i 



-M T 



M — , 

2kq T Jo 



dt x ■ ( 



(55) 



r=o 



As we have seen in Eq. ([38]) the supersymmetry of the action 5* leads to the result that the 
expectation value of x ■ ( is time-independent and thus can be evaluated at any time t, in 
particular at t — T. We then can perform the t-integral and obtain 



G'(x,y) 



dT N(T) exp 



2k q 



-M T 



2k 



M-ix(T)-((T) 



(56) 



where the average with respect to the action S is defined in Eq. (^). For the calculation of 
the last average we use the well-known fact (see e.g. Ref. J28|) that the expectation value 
of time-ordered products of Heisenberg operators Onit) = exp(iHt) O exp(-iTCt) is given 
by the insertion of 0(t) in the corresponding path integral. Thus 



< x,T\ T 



x H {h 



\y, >= < x\ e 



-iHT 



T 



\y>= (x{t 1 )-C{t 2 )) s 
(57) 



since Eq. ( |17|) tells us that the (Weyl ordered) 7 matrices are to be replaced by 2£. Differ- 
entiating with respect to ti (the equal time contribution vanishes) and putting t\ = t 2 = T 
we obtain 



<x(T)-({T) > s = <x\i H,x, 



-L e~ mT \y > . 
2 |y 



(58) 



Evaluating the commutator with the help of Eq. flTOf) and the canonical commutation rela- 
tions we find 



< x(T) ■ C(T) > 5 



1 



2k 



< x\ f[ e - lH1 \y>= — )p x + ig4(x)\ < x\ e~ lH1 \y > (59) 



which, inserted into Eq. (|56|) , gives exactly the same result for the propagator as the global 
projection method, i.e. 



G'(x,y) = G(x,y). 



(60) 
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IV. NONRELATIVISTIC LIMIT 



If the mass M of the fermion becomes large the integral over the proper time T is 
dominated by the stationary points of its integrand which approximately occur at IIq = M 2 . 
Therefore we make the Ansatz 



n = sM + E , s = ±1 . 
For the nonrelativistic limit it is very convenient and natural to take 

kq = M 

and to assume 

e = o : 1 

M 



(61) 
(62) 

(63) 



In this Section we will write V d x, V d p and iV d spm with d = 3, 4 to stress the different 
dimensionality of relativistic and nonrelativistic path integrals. In the global projection 
method we then obtain from Eq. fli~7|) 



G(x,y) 



2M 



^-^ + M)expf T -— j ]T / °° rfTexp (-i MT ) 
• iV 4 spin J V 4 x V 3 U V% J VE exp{ tS s [x, U, E, Q } r=Q 



(64) 



with 



s s [x,n,E,(] = - i ao)-aT) + 



dt 



M 



— Msxo + II • x — gA xo 



+gA-k 



ig 



+ iC-C- —F UU CC + sE- Ex + 

According to our assumption the last term in the square bracket is 0(1/ M 3 ) which we 
neglect. The path integral over E then gives a functional 5-function 



E 2 



(65) 



N 



S \xq — s] — lim TT 5 . 



/ x 0,k — %Q,k-l 



At 



At 



N 



(66) 



The functional integration over x can now be performed trivially, with the result that the 
time co-ordinate has the proper time dependence 



xo(t) = yo + st . 



(67) 



However, one 5-function remains because there are only (N—l) integrations in the discretized 
path integral for the co-ordinates: 



G(x,y) 



2M 



g4 + M)exp^y~ y j £ J™ dT 5 (x -y - sT) exp (-iMT) 



■ Nl pin 



V 3 xV 3 UV\ exp C(0)-C(T) 



dt 



i( ■ c + n • x 



2M 



gsA + gA-±- ^F^C 



(68) 



r=o 
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The remaining 5-function enforces the boundary condition xq(T) = xq and can be used to 
perform the integration over the proper time T yielding 



T 



s {xq - y ) 



(69) 



In other words, in the non-relativistic limit the proper time becomes the ordinary time (dif- 
ference), as expected. Since the proper time is positive, the (s = +l)-term describes forward 
propagation of the particle whereas the (s = — l)-term describes backward propagation of 
the antiparticle, which is also contained in the Feynman propagator but decouples in the 
heavy mass limit. Furthermore, the global projection operator in front of the propagator 
(ESI) can be replaced by 



2M 



{0 X -g4 + M) 



1 

M 



(70) 



as the xo-derivative acting on the phase factor exp (—iMs(xo — yo)) gives the leading con- 
tribution. Since 



To 



1 
-1 



(71) 



the (anti-)particle propagator acts only on the (lower) upper components of Dirac spinors if 
the remaining path integral is diagonal in 2 x 2 Dirac space (which will turn out to be the 
case). Shifting back to integration over p we therefore obtain 



G(x,y) 



9(s(x -y )) -(l + s7o)e 



-iMs(x -yo) 



exp 7 • 



s=±l 



or 



N Bpm 



I 



T=s{xa-yo) 



V 3 x V 3 p V\ exp | C(0) -C(T)+i 
1 



dt 



*C • C + P - * 



2M , P -9Af- 3 sA -^ CC 



r=o 



The time- dependence of the electromagnetic potentials and fields is fixed by Eq. 
stituting 

t' = y + st, t' E [yo,x ] 
x(t) = x'(t') , p(t) = p'(t') , C(*) = C(t') 



(72) 
Sub- 

(73) 



the boundary conditions for the co-ordinate space path integral become the usual ones for 
a nonrelativistic path integral |29|] 

(74) 



x'(z/ ) = y> x '(^o) = x. 

Omitting the primes, the action in the phase space path integral now reads 



S s [x,p,C] = -iC M (2/o) C(^o) + 
(P - ^A) 2 



dt 



^Cu + p-x-W s (x,p,2C) 



K(x,p, 2C) = 8 



M + 



2M 



+ gA + l j^F„ u CC ■ 



(75) 
(76) 
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Here we have absorbed the phase factor exp(— iMs(xo — yo)) into the Hamiltonian 7i s . 
Finally we simplify the spin degrees of freedom by using 

F^C = 2Co E • C - B • (C x C) (77) 

and observing that the first term in Eq. (|77f) is linear in ( . After shifting ( Q = T /2 + £ , 
the ^-integration can be performed and leads to a term in the remaining action which is of 
0(1/M 2 ) : the Fourier transform of a Gaussian is again a Gaussian. In that process part of 
the spin normalization factor is cancelled. The same argument can be applied with respect 
to the £ - integration so that in leading order only igsToE ■ T/M survives from the first 
term in Eq. (|77|). Performing the required differentiations with respect to To we see that 
we obtain a contribution to Eq. ([?(]) of the same order which was already neglected. Notice 
that 2£ E • £ is "odd" in the sense of connecting large and small components in the Dirac 
equation; from the standard Foldy-Wouthuysen transformation of the Dirac Hamiltonian it 
also follows that the odd parts are suppressed by a factor 1/M compared to the "even" ones. 
In addition, since 








-<7i 




(78) 



one can set 



exp (7 • V r ) 



exp [icr • Vr) 



(79) 



for the remaining even part of the action. After changing the signs of the spin terms by the 
substitution £ — > iC, the nonrelativistic limit for the Dirac propagator finally becomes 



G(x, x ;y,y ) 



with 



s=±l 
x(x )=x 



~ £ 6 (s(x - y )) - (1 + s l0 ) e CT - Vr 



j V 3 xV 3 p J V 3 ( exp|c(z/o)-C(xo) 



*(yo)=y 



+i / dt 

Jyo 



C(s/o)+Cor )=r 

zC-C + p-x-H s (x,p,2C) 



r=o 



(80) 



N spm 



V 3 C exp 



C(Vo) ■ CM 



.C(w)+C(^o)=r 



The Hamiltonian 

K(x,p,2C) = s 



M + 



(P - <?A) 5 



2M 



X() 



dtC-C 



wo 



(81) 



+ ^ + |B-(CxC) 



(82) 



coincides exactly with the standard Foldy-Wouthuysen Hamiltonian for particles and an- 
tiparticles (see, e.g. Ref. Hl9| , Eq. (4.5) ) 
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7o 



M + 



(P ~ flA)' 
2M 



+ flAi — 7n CT ■ B 



if the last term (the so-called Pauli term) is rewritten using 

(T x it = 2z er . 



(83) 



(84) 



This ensures that the time evolution is governed by a Grassmann-even Hamiltonian. It 
is, of course, straightforward to start from the nonrelativistic spin-dependent Hamiltonian 
(jB3|) and by using the identity fl84|) to derive the 3-dimensional path integral (|80| ) for the 
propagator as has been done in Ref. |J, Section 5. Here we proceeded in the reverse order 
showing how the nonrelativistic limit can be taken within the path integral representation 
of the Dirac propagator. It should also be possible to evaluate higher-order terms in the 
nonrelativistic reduction in this way or to obtain the semiclassical limit of the propagator 



V. DIRAC PROPAGATOR IN AN EXTERNAL SCALAR POTENTIAL 

For some applications one needs the propagator of a fermion which moves in an external 



scalar field S(x) as well. For example, in the Walecka model |3T| the exchange of a scalar 
meson generates attraction between nucleons whereas massive vector mesons are responsible 
for repulsion at shorter distances. In such cases we need to evaluate the following Green 
function 



G(x,y) 



where 



1 



H- M*(x) +i0 



(85) 



M*(z) = M + S(x) 

is the effective, position-dependent mass of the fermion. 

The previous method of multiplying numerator and denominator in Eq. 
obviously does not work anymore since 

(H- AT) (H+ M*) = H 2 - M* 2 + [H, M*] 



(86) 



by H+ M* 



(87) 



is not Grassmann even. Consequently there are statements in the literature [27] that in this 
case a five- dimensional formalism is the only possible approach. However, this is not the 
case: the problem to rationalize the denominator of the Green function is analogous to the 
problem of inverting complex matrices by using only real arithmetic. This is easily achieved 
by writing 



1 



A + iB 



(l-iA^B 



1 



A + BA- X B 



(88) 



Therefore we have 
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p.- M* V M + M*^) flj^l/L- MM* + iO ^ 

as the Di Vecchia-Ravndal representation for the present case. Again we have the choice to 
project on fl= M* either during the evolution or at the end. If we adopt the latter approach 
the quantum mechanical Hamiltonian which governs the proper time evolution is now 

W = — ft,, (9°) 

2/V M*(xy v ; 

The phase-space path integral representation of the propagator is now determined by 
evaluating the Wigner transform of Eq. (|90| ) (see Appendix A). With the abbreviation 
U(x) = M/M*(x) one obtains 

H{p,x,i) = -^-Tl 2 U(x) + ^ lfllu F^(x)U(x) 
2kq Akq 

-T-Wv [d»U{x) IT - Wd"U(x) ] - J-d 2 U(x) . (91) 

Since the Hamiltonian is quadratic in II = p — gA the momentum path integral can still be 
performed so that the Lagrangian path integral representation of the propagator reads 

iMM*(x)T\ ( d\ 



G(x,y)=-^-[U(x)(0-g4(x)) + M} jT dT N(T) exp ^ j exp (7- 

px(T)=x /• ( pT 

/ Vx V( exp C(0) • C(T) + i / dt L(x, x, (, () (92) 

Jx(0)=y J JO I n f\ 



C(o)+c(T)=r 

with 



r=o 



K .2 4/ x ■ *0 

1 „„, , 1 



L(x, x, C, C) = «C • C - 7777^ - ^ (x) • x - ^f/(x)F^(x)C M C 

2(7 (a;) 



-2ix ■ CTrfrC • <9t/(x) + —d 2 U(x) + 2i5(0) log U(x) . (93) 
(7 (x) 8^o 



The last term arises from the quadratic fluctuations 



n^T^y = exp^-2]Tlogf/(:r fc )) = exp (i2i ± At^logC/(x^ 



(94) 



in the discretized momentum path integral which are now position-dependent due to the 
effective mass M*(x). The awkward 5(0) appears as the formal limit of 1/At when the 
time-slicing At is made infinitesimal (see Ref. [32j| , Chapter 19) and cancels consistently 
against other divergencies 



VI. EFFECTIVE ACTION FOR QUENCHED QED 

In order to reduce the number of degrees of freedom it is advantageous for some appli- 
cations to integrate out the bosons which mediate the interactions. The price to be paid 
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is, of course, a more complicated two-time effective interaction. We will outline this pro- 
cedure by considering Quantum Electrodynamics (QED) (or the Walecka model without 
scalar mesons ff) 



C = £o(A)+$(i9-g4-M )il> 



(95) 



where 



C (A) = -if^ + l -m 2 A 2 -~A(d • A) 2 



(96) 



is the Stiickelberg Lagrangian with a gauge parameter A. We have given the photons a mass 
m in order to regularize infrared divergencies. 

The generating functional for the 2-point function with an arbitrary number of photons 



is 



Z' \j, x] 



VA ( x 



i@-g4- M 



0) exp{tA [A} + (j,A)} . 



(97) 



Here the free vector meson action is denoted by *4.o[^4] = / d A x Cq(A) and we have neglected 
closed fermion loops (quenched approximation) j34j in order to have a single world-line for 
the fermion. For integrating out the vector field A, we use the path integral representation of 
the Dirac propagator in an external vector field in its supersymmetric form and the gaussian 
integration formula 



J VA, exp [i (A„ {G~ l r v A v ) + (Ap, h») 



oc exp 



Here 



GT{k) 



gf iu _ j m 2 ^ k^k v /m 2 



k 2 — m 2 



k 2 — m 2 /X 



(98) 



(99) 



is the standard propagator for massive vector particles (Ref. @]). From the linear terms in 
A^ we read off 



K{y) = JM +ig dt d6 DX^t, 6) 8\ y - X(t, 9) ) . 



(100) 



For the present purposes the local projection method is preferable because the whole 
dependence on the photon field resides in the free photon action and the electron-photon 
interaction. After integration over the photon field we then obtain for the generating func- 
tional (191) 



4 If one also wants to integrate out the scalar mesons, the 5-dimensional Berezin-Marinov descrip- 
tion has to be used because only then a Gaussian path integral for the scalar mesons is obtained; 
our 4-dimensional form (^) is highly nonlinear in S(x) . 
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Z' [j, x] = const, exp 7 



' dT J Jo 

dx exp 



dT N(T) exp 



2k 



-M 2 T 



2k q 



Mx VxV( e 



iS eS [X,j] 



(101) 



r=o 



where the effective action is given by 



S eS [X,j) = S '[X] + - (hp, G» v h h 



(102) 



As in Ref. [j34j| , it is advantageous to split it up into terms involving zero, one or two external 
sources j(y)- The latter one leads to disconnected diagrams and can be discarded. We then 
have 



S eS [X,j] = S '[X] + + S 2 [X,j] 

where the free action is given in Eq. (H3) and the interaction part by 



(103) 



Si[X] = ^ [ dt x J dB x dt 2 J d6 2 J -^-G^(k) DX^ti, 9i) DX v (t 2 , 9 



d A k 



2 Jo 

Note that the "current" 

J,{X) = DX^t,. 



exp{-ik-[x(t 1 ,e 1 )-x(t 2 ,e 2 )]} . (104) 



-6xJt)+a(Jt) 



(105) 



looks like in scalar QED but is Grassmann odd and does not depend on the integration 
variable k. Therefore the ^-integration can be performed easily giving the photon propagator 
in configuration space with argument X(ti, 61) — X(t 2 , 9 2 ) Written in components the 
interaction term 



2 



d A k 



4 cr{k) 



i M (ti) + — Ut^k-Cih) 

K 



Xv(t 2 ) (v(t 2 ) k ■ C(t 2 ) 

K 



-ik-[x(ti)— xfo) ] 



(106) 



is seen to contain up to quartic terms in the spin variable (. This means that, unlike the 
case of external fields, the Grassmann variables cannot be integrated out anymore to give 
a "spin factor". Vice versa, it is impossible to eliminate the photon field starting from the 
spin factor formulation for the propagator. 
The source term becomes 



S 2 [X,j] = ig j cfiyjM J q dt jd9 j — G^(k) DX v (t, 9) exp { -ik ■ [X(t, 9)-y]}. 

(107) 



It is also possible to use the global projection method which does not have a spin-orbit 
coupling. However, there is an additional dependence on the photon field in the covariant 
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derivative acting on the path integral in Eq. fll8|) which makes it less suitable for deriving 
an effective action. 

To conclude this Section, we note that the effective action in Eq. ( |f03| ) allows a particu- 
larly concise derivation of the transformation properties of Green functions f\ under a change 
of the gauge parameter A : We see from the photon propagator G^ u {k) in Eq. ([Mj) that 
a change in A only effects the term proportional to k fl k u . For this term the integrals over 
the proper times ti and 6i occurring in the effective action may be performed exactly as the 
integrand is a total derivative, i.e. 



dt ddk- DX(t, 9) e 



-ik-X{tfi) 



dt dOD e 



-ik-X(t,6) 



(1 



— e 



-ik-x 



(108) 



The change in Si (Eq. ( |104| )) and S 2 (Eq. ( |107| )) induced by a change in A from Ai to A 2 , 
say, is therefore only dependent on the endpoints of the path x(t) and not on the path itself. 
If we define A(x 2 ) to be the Fourier transform of the change of the coefficient (= A(k 2 )) of 
kpky in the photon propagator, i.e. 

d 4 k ( 1 1 N 



A(x 2 



77T 



-ik-x 



(2vr) 4 \k 2 -m 2 /\ 2 ' k 2 -m 2 /\ y 



(109) 



then the corresponding change in S\ is given by 

g 2 f d 4 k 
~2 



6S 1 = S( 



A=A 2 



[X] - S f*[X\ = 9 - A(*V (W**) (l-e* 8 



9 



A(0) 



while the change in is 

5S 2 = ~ig 



d 4 yj(y)-d y [A([y 



A(y 2 



(110) 



111) 



Note that not only is the change in Si^ independent of the path, so that it may be pulled 
out of the path integral in Eq. ( |101| ), it also does not involve the Grassmann valued T nor is 
it dependent on the proper time T. Hence the generating function for the Green functions 
with gauge parameter \ 2 is related to that with gauge parameter Ai in a very simple way, 
namely 

Z' X2 [j,x]=e^ +ss ^Z' Xi [j,x] (112) 



exp 



{ ig 2 [A(x 2 ) - A(0) ]+gJ d*yj(y) ■ d y [ A ([y - x} 2 ) - A(y 2 ) } } Z' Xi [j,x}. 

As special cases we can derive the transformation laws for the propagator and the electron- 
photon vertex from this expression Setting j = we obtain 



5 These transformations were first derived for the electron propagator and the electron-photon 



vertex by Landau and Khalatnikov [35] and extended to general Green functions by Fradkin and 
Zumino |36f. 



6 See Ref. [35]. Note that in that paper the photon propagator is defined with a minus sign with 
respect to ours. Hence our function A(y 2 ) is —Ap(y) of Ref. [35] and our untruncated vertex 



function is the negative of the function Bn defined by Landau et al. 
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;ii3) 



while by differentiating once with respect to the current and then setting j = we obtain 
the (untruncated) vertex function 



(114) 



It should be noted that these relations are valid even if the photon mass (which violates 
gauge invariance) is kept non-zero in the photon propagator. 



VII. SUMMARY AND CONCLUSIONS 

The main purpose of this work is to explore a four-dimensional path integral representa- 
tion for the Dirac propagator general enough to describe the particle's motion in both vector 
and scalar fields. Although the four- dimensional approach, for an external vector field, was 
proposed by Di Vecchia and Ravndal almost twenty years ago, it had received limited at- 
tention up to now. Instead it is standard to use the Berezin - Marinov approach where one 
introduces a fifth component to eliminate the extra spin degree of freedom. However, the 
fifth component has no clear physical meaning and the necessity of introducing a 75 for the 
evaluation of the propagator seems rather unnatural. The four-dimensional representation 
avoids these difficulties; in addition the supersymmetry transformations become easier and 
more natural to generate. 

Working within this four-dimensional formalism we have presented two alternative meth- 
ods to project out the unwanted spin degree of freedom. The first method projects onto 
the final state after the time evolution and is hence termed global, whereas in the second 
method the projection is done at each step in the time evolution and it is therefore referred 
to as local. Extending previous work by Reuter, Schmidt and Schubert we have shown that 
due to the supersymmetry the two methods are completely equivalent and may be used 
according to convenience. The main difference between the two approaches is that the path 
integral representation using the local projection has an explicit spin-orbit coupling term. 
It was therefore crucial for the proof of equivalence to generalize the results of Ravndal and 
Di Vecchia regarding the supersymmetry transformations to apply also in the case where 
the spin-orbit term appears. In Refs. [0,|T31 it was pointed out that, in the case where no 



spin-orbit term was present, a simple way of generating both bosonic and fermionic trans- 
formations is to shift the times t and 9. We show in this paper that in the presence of a 
spin-orbit term in addition to a shift an appropriate scaling of the times t and 9 is needed 
in order to generate the correct supersymmetry transformations. This scaling is such that 
the parameter k and the supertime x remain unchanged. 

For the case of a Dirac particle in an external scalar potential it was generally believed 
that a five-dimensional approach was unavoidable. We have here shown that this is not 
the case and we used the four-dimensional description to obtain the Dirac propagator in an 
external scalar field . 

Despite the attention given to spin in the path integrals, a nonrelativistic reduction 
starting directly from the Dirac propagator was still missing. By expanding the relativistic 
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expression in powers of 1/M we were able to reduce the path integrals to three-dimensional 
form and to obtain the leading nonrelativistic result described by the Foldy-Wouthuysen 
Hamiltonian. 

Finally we applied the four- dimensional approach to quenched QED in order to obtain 
a supersymmetric formulation for the generating functional of Green functions with one 
electron line and an arbitrary number of external photon lines. It was possible to do this as 
in the quenched approximation the photons can be integrated out, yielding a path integral 
only in the electron degrees of freedom, albeit with a complicated nonlocal interaction. In 
this form one can apply methods along the same lines as those used in the study of the 



polaron problem as described in Ref. |2q| . Furthermore, we showed that it is a rather 
simple matter to derive the Landau-Khalatnikov transformations for the propagator, vertex 
function and indeed any higher-point function from this formalism. 

From a field-theoretic point of view, the worldline technique is particularly appropriate 
whenever one deals with a situation where internal fermion loops may either be neglected 
or taken into account perturbatively. As this situation arises quite naturally in the non- 
relativistic regime, the technique would appear to be particularly appropriate in that setting. 
We think that the reason it has not received a great deal of attention by physicists working 
in that area is partly due to the fact that the commonly used five- dimensional representation 
appears artificial within this context. In this paper we have tried to convey the message that 
for most problems the five-dimensional formulation is not only unnecessary but in fact less 
transparent than the four- dimensional one. It is our hope, therefore, that this paper makes 
worldline techniques more accessible to a wider audience than they have been up to now. 
Note added in Proof: Recently it was pointed out to us that the elimination of the fifth spin 
variable was also considered by T. Allen using Hamiltonian methods (T. Allen, Phys. Lett. 
B 214, 87 (1988); see also T. Allen, Caltech Ph.D. Thesis, (1988) ). Also, J. W. van Holten 
has advocated the use of a commuting rather than anti-commuting fifth spin variable and a 
different four-dimensional approach (see the first reference in [10] as well as a more concise 
discussion of the problem in Nuc. Phys. B (Proc. Suppl.) 49, 319 (1996) ). Finally, another 
important contribution to the literature on spin in path integrals missing from Ref. [6] is 
the paper by M. Halpern, A. Jevicki and P. Senjanovic, (Phys. Rev. D 16, 2476 (1977) ). 
We are grateful to Profs. T. Allen, M. Halpern and J. W. van Holten for correpondence 
regarding these references. 
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APPENDIX: SPIN PATH INTEGRAL FOR THE TIME EVOLUTION 

OPERATOR 

Here we consider the matrix element of the time evolution operator 

U(x,y) = (x I exp {-iKT\ I y) (Al) 
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where 



(A2) 



is a Weyl-ordered Hamiltonian. Breaking up the evolution operator in N time steps we 
obtain in the usual way 



U(x,y) = lim / d 4 Xi...d A x^_i 

N^oo J 



d A pi d A p N 



(2tt) 4 (2tt) 4 
exp [ -iTCwipN, x Ni Tat) At ] 



exp 



N 

i ^ t Pi ' i^i ^d—l, 
i=l 

exp [ -iHwipi, x 1 ,ji)At 



with xq — y and xn = x. Here 



d*y (x — ^- 



n 



X 2 



ip-y 



(A3) 



(A4) 



is the Wigner transform (or Weyl symbol) of the Hamiltonian which is the closest classical 
analogue to the (Weyl-ordered) quantum operator P7[| . We will suppress the subscript "W" 
in the following. 

There are two essential steps to derive a path integral with spin: 

(i) Because the Dirac matrices do not commute the ordering of the factors is essential and 
the exponentials cannot be combined with impunity. As is well known this also happens 
in ordinary quantum mechanics for time-dependent Hamiltonians. We therefore have 
assigned an artificial time-dependence to the Dirac matrices and can write now the 



time evolution operator as a time-ordered path integral [29 

U(x,y) = J VxVp T exp j-i £ dt[p ■ x + H(p(t),x(t),j(t))] 

5 \ 



VxVp exp 



r-T 

ii dt 
h 



p- x + H\ p(t),x(t), 



hit), 



Texp 



dtff{t)^{t) 



(A5) 



J pf=0 



Here p^{t) are Grassmann sources which are assumed to anticommute with the Dirac 
matrices. The boundary conditions for the x-space path integral are 



2^(0) = y^ 



x^T) 



x 



(A6) 



The time-ordering symbol T would be disastrous for further manipulation of the path 
integral. However, in the special case it can be eliminated by the relation 



V(T) = Texpl£ dt 



exp 



T 



dt! / dt 2 ff{t x )p^t 2 ) 



•<; / n 'I' i>"{>)' 



i-i- 



(A7) 
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This can be proved by solving the corresponding evolution equation 



dV(T) 
dT 



p»{T) llx V(T) , V(0) 



(A8) 



by using the Magnus expansion [33 



V(T) = expjj^ dt(f{t) llx + \ f dU f 1 dt 2 [p M (ti)7 M ,P,(t 2 )7l +••• 



2 Jo 



• (A9) 



The commutator yields —2p fl (ti)p ll (t2) which is a commuting c-number so that all 
higher terms in the expansion which involve multiple commutators vanish . On the 
right-hand-side of Eq. ( |A7|) we can now drop the artificial time- dependence of the 
Dirac matrices. 



(ii) The differentiations with respect to //*(£) which are required in Eq. (^5) can only be 
performed easily if they appear linearly in the exponent. This can be achieved by 



"undoing the square", which is a standard procedure f29 |. However, because p^(t) is 
anticommuting and one needs an even object in the exponent as evolution operator, 
we have to do it with the help of a Grassmann path integral. We thus use the identity 



expj-^ dti jT dt2p t *{ti)p fl {t 2 ) 



T 



Vi exp / dt -Ut)e(t) + 2f^(t)Ut) 



. (A10) 



and the antiperiodic boundary condition £^(0) + ^(T) = for the Grassmann path 
integral. The standard way of proving this identity in the continuum formulation is 
by solving the (differential) equations of motion which should give the exact result for 
quadratic actions. However, it is very useful (and reassuring) to have an unambiguous 
formulation with finite time steps At, which we will present now: the discretized form 
of 



S = - f dt \-^(t)e{t) + 2(f{t)Ut) 

J 



(All) 



may be written as 

N 



s 



E 

i=i 



(A12) 



where At = T/N and N needs to be even for the path integral to be an even quantity. 
In discretized form the path integral over £ in Eq. (|A10 ) can now be done by the 
stationary phase method. The (difference) equation of motion 



fc+i 



pk-i 



(A13) 
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can be solved using antiperiodic boundary conditions for £, i.e. £^ = — £o > £jv+i 



It is convenient (but not necessary) to impose the equivalent boundary conditions for 
p. Note that the particular discretization of p(t) ■ £(£) in Eq. ( |A12| ) is chosen so that 
the equations of motion (|A13|) for the odd and even sites are coupled. This avoids the 
infamous "fermion doubling" problem. The solution to the equation of motion is 



At j At N 

z fe=i z k=i 



(A14) 



Substituting the solution £ cl into Eq. ( |A12j ) we find 



N 



i-1 



i=l k=l ° i=l 

The path integral over £ can now be performed yielding 

/P£exp[-S] 



expf-Sc 



(A15) 



(A16) 



since the determinant from the quantum fluctuations is canceled by the denominator. 
Taking the continuous limit of S c \ only the first term in Eq. (|A15|) survives and we 
obtain the required result. Having proven the relation (|A10|) by writing the functional 
integrals in a well defined discretized form we can now use it in Eq. ( |A5| ) with all 
manipulations formally done in the continuum. 



Using the representation 

exp| [ T dt p»{t) lfl 



o 



(A17) 



r u =o 



we obtain 



U {x, y) = exp 7 



Here 



or 



VxVpV£N spm exp{-i 



dt 



p-x-i£-£ + H(p,x 7 2£ + r) 



N 



spin 



n -1 



T>£ exp 



dtLe 



r=o 
(A18) 



(A19) 



is a normalization factor for the spin integral. Note that the operation in Eq. (|A17|) is 
in general not just a replacement of the boundary variable T by the corresponding Dirac 7 

7/^ ; 



matrix but involves an antisymmetrization as well. For example, exp {7 ■ d/dT} 
but 

2 



r=o 



exp < Y 



d 



r r 



r=o 



7 • 



d_ 
dT 



r r 



r=o 



7 • 



d_ 

dr 



7^ + T^ v 



( + liflv ) 



r=o 
(A20) 



This is the inverse transformation of the Weyl representation for fermionic operators. 
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